Partial separability revisited: Necessary and sufficient criteria 
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We extend the classification of mixed states of quantum systems composed of arbitrary number of 
subsystems of arbitrary dimensions. This extended classification is complete in the sense of partial 
separability and gives 1 + 18 + 1 partial separability classes in the tripartite case contrary to a 
former 1 + 8 + 1. Then we give necessary and sufficient criteria for these classes, which make it 
possible to determine to which class a mixed state belongs. These criteria are given by convex roof 
extensions of functions defined on pure states. In the special case of three-qubit systems, we define a 
different set of such functions with the help of the Freudenthal triple system approach of three-qubit 
entanglement. 
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I. INTRODUCTION 

The notion of entanglement [l|, was regarded by 
Schrodinger 0, 0| to be the characteristic trait of quan- 
tum mechanics. It serves as a resource for Quantum In- 
formation Theory a relatively new field of research 
dealing with the properties, characterization and appli- 
cations (mostly in quantum computation [5] ) of the non- 
local behavior of entangled quantum states. 

For a multipartite quantum system being in a pure 
state, it is easy to decide, in general, which subsystems 
are entangled with some of the others or, equivalently, 
which subsystems can be separated from the others. For 
a multipartite quantum system being in a mixed state, 
however, this partial separability problem has not been 
considered in the full detail yet. This problem is twofold. 
Even if we have the definitions of the different classes, 
which is not self-evident at all for more-than-two-partite 
systems, deciding to which class a given state belongs is 
also a nontrivial task. In this paper, we work out solu- 
tions for both parts of this problem. Considering the first 
part of this problem, we extend the classification based 
on /c-separability and afc-separability given by Seevinck 
and Uffink [6( , which is the extension of the classification 
dealing only with afc-separability given by Dur and Cirac 
BUI- We discuss in detail the tripartite case, then give 
the definitions for systems of arbitrary number of subsys- 
tems. Before we outline our solution for the second part 
of the problem, we take a short detour. 

If in the tripartite case we restrict ourselves to qubits, 
which is a relatively well-understood chapter of the the- 
ory of quantum entanglement, some interesting results 
are known from the literature. From the point of view 
of the present work, the most important ones are the 
following three. First, (i) this is the system where the 
nontrivial structure of entanglement manifested itself for 



the first time and it came to light that "there are differ- 
ent kinds of entanglement" of pure states [9|. Then, (ii) 
these different kinds of pure-state entanglement Q give 
rise to classes of mixed state entanglement [lfj. On the 
other hand, (iii) recently a beautiful correspondence was 
found between the three-qubit Hilbert space and a par- 
ticular FTS (Freudenthal Triple System), a correspon- 
dence which is "compatible" with the entanglement of 
pure three-qubit states [llj. Apart from these three, for 
the sake of completeness, we have to make mention of (iv) 
the famous phenomenon of monogamy of qubit systems, 
which was revealed first for three-qubit systems [l2| then 
shown for multiqubit systems [13| . and (v) the interest- 
ing twistor-geometric approach of the entanglement of 
three-qubit systems [14] . 
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Item (iii) above gives us a hint of an answer to the sec- 
ond part of the problem for three-qubit systems. In the 
FTS approach of three-qubit entanglement some special 
quantities have appeared, from which we gain real val- 
ued functions on pure states. These functions have very 
useful vanishin g prop erties, which enable their convex 
roof extensions [15|-[l7| to identify all the classes that our 
extended classification deals with. On the other hand, 
it will be possible to define suitable functions for the 
identification of the classes in general, for subsystems 
of arbitrary dimensions — moreover, for arbitrary num- 
ber of subsystems — in another way than was done with 
the FTS approach working only for three qubits. How- 
ever, we will keep the considerations coming from the 
FTS approach, because these considerations have given 
us the main ideas, they have advantages for the case of 
three-qubits, and, besides these, they are beautiful and 
interesting in themselves. 

In the bipartite state — either pure or mixed — 

can be either separable or entangled [l8[ , and the vanish- 
ing of the convex roof extension of local entropies of pure 
states is a necessary and sufficient criterion of separabil- 
ity. For us, this is the archetype of the general method 
of the detection of convex subsets by convex roof exten- 
sions. However, for more-than-two-partite systems, the 
partial separability properties have a complicated struc- 
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ture, and, to our knowledge, this method was not used. 
Instead of that, the usual approach was the use of witness 
operators, as was done originally for three-qubit systems 
[10|, or other necessary but not sufficient criteria for the 
detection of convex subsets [jj, |2], LUJ |20| . 

Before starting, we review the classification schemes 
of states of multipartite quantum systems. One of the 
main concepts here was the use of LOCC (Local Oper- 
ations assisted by Classical Communication [2l|) either 
with certainty or with possibility, for the purpose of clas- 
sification. This concept has turned out to be useful in 
the restricted case when the input and output states are 
both pure. First we recall the classification schemes deal- 
ing with LOCC. For mixed states, only coarse-grained 
classifications are worked out, which are recalled as well. 

LOCC classification: Two states are equivalent under 
LOCC — they are in the same LOCC class — by definition 
if they can be transformed to each other with certainty 
by the use of LOCC. For pure states, it turned out that 
two states are equivalent under LOCC if and only if they 
can be transformed into each other by LU (Local Uni- 
tary) transformations [21] . So, for pure states, this gives 
the most fine-grained classification scheme imaginable. 
Many continuous and discrete parameters are required 
to label the LOCC classes |22H26j . From the point of 
view of quantum computational purposes, two LOCC- 
equivalent pure states can be used for exactly the same 
task. However, to our knowledge, there is no such prac- 
tical criterion of LOCC equivalence for mixed states as 
the LU equivalence was for pure states. 

SLOCC classification: A coarse-grained classification 
can be defined if we demand only the possibility of 
the transformation. Two states are equivalent un- 
der SLOCC (Stochastic LOCC) — they are in the same 
SLOCC class — by definition if they can be transformed 
into each other with non-zero probability by the use of 
LOCC. For pure states, it turned out that two states are 
equivalent under SLOCC if and only if they can be trans- 
formed into each other by LGL (Local General Linear) 
transformations @. (Sometimes that was called ILO, 
stands for Invertible Local Operation [gj], but we prefer 
the uniform naming after the corresponding Lie groups.) 
So this gives a coarse-grained classification scheme for 
pure states. In some cases, including the three-qubit case, 
only countable finite SLOCC classes arise [9(. From the 
point of view of quantum computational purposes, two 
SLOCC-equivalent pure states can be used for the same 
task but with a different probability of success. Again, 
to our knowledge, there is no such practical criterion of 
SLOCC equivalence for mixed states as the LGL equiva- 
lence was for pure states. 

PS classification (Partial Separability): A more 
coarse-grained classification involves only the partial- 
separability properties. This works for both pure and 
mixed states and gives only countably finite classes in 
both cases. We elaborate this classification in detail in 
this paper for mixed states. This classification deals with 
all the possible kinds of partial separability, which are of 



finite number, whose special cases are the subsets of k- 
separability and o^-separability @. From the point of 
view of quantum computational purposes, however, this 
classification is a bit too coarse grained, since it does not 
make distinction among pure states contained in differ- 
ent SLOCC classes but having the same PS properties, 
although these states may be suitable for different tasks. 

PSS classification (Partial Separability extended by 
pure-state SLOCC classes): A cure for the problem above 
is another means of classification, which was given by 
Aci'n et. al. [l(| only for three-qubit states. Here, the 
starting point is the pure-state SLOCC classes which are 
of finite number, and the only difference between the par- 
tial separability classes and SLOCC classes is the split of 
the three-qubit entangled class into two classes [9(. The 
PSS classes arising from these classes for mixed states are 
the same for biseparability, and only the tripartite entan- 
gled set is divided into two classes. This classification has 
the advantage of differentiating among different SLOCC 
classes of pure states, and also among mixed states de- 
pending on which kind of pure entanglement is needed for 
the preparation of the state. However, in the majority of 
the cases there are continuously infinite SLOCC classes 
of pure states labeled by more than one continuous pa- 
rameter [9|, l27|, |28[ , in which case it is not clear how this 
classification can be carried out, if it can be at all. 

The organization of this paper is as follows. In the 
first half of the paper, we work out the main concepts 
on three-qubit states. In Sec. [H] we review the SLOCC 
classification of pure three-qubit states. We recall the 
conventional LU invariants (in Sec. Ill Ap and the LSL 
tensors (Local Special Linear) of the FTS approach (in 
Sec. |HB|) by which the SLOCC classes can be identified. 
Then we obtain a new set of LU invariants (in Sec. Ill C[) 
being necessary later for mixed states. In Sec. IIIII we 
elaborate the PSS classification (which contains also the 
PS classification) for mixed three-qubit states. We de- 
fine the PS(S ) subsets (in Sec. |HT~A"|) and PS(S) classes 
(in Sec. IIIIBI) . Then we give the functions for the iden- 
tification of the PS(S) classes (in Sec. ITHCl In Sec. [TV] 
we demonstrate the nonemptiness of some of the new 
classes for the three-qubit case by explicit examples. In 
Sec. El we generalize the functions for the case of three 
subsystems of arbitrary dimensions. First we see how 
far the method coming from the FTS approach can go 
(in Sec. IV Al) ; then we formulate a more general set of 
functions working without limitations (in Sec. IVBp . In 
Sec. I VII we generalize the construction for the case of 
arbitrary number of subsystems of arbitrary dimensions. 
We work out the labeling of the PS subsets (in Sec. I VI A)) 
along with the PS classes and give a general conjecture 
about their nonemptiness (in Sec. IVIBI) . Then we con- 
struct the functions identifying the PS subsets and classes 
with the minimal requirements (in Sec. IVI Cp . as well 
as with stronger requirements leading to entanglement- 
monotone functions (in Sec. IVIDj) . In Sec. IVIII we give 
a summary, some remarks, and open questions. Some 
technicalities about the new set of three-qubit LU invari- 



ants and proofs of some statements about the general 
construction are left to Appendixes [A] and [B] 



II. PURE THREE-QUBIT STATES 

Before starting, we set some conventions that are very 
convenient for the tripartite case. The labels of the sub- 
systems are the numbers 1, 2, and 3, while the letters 
a, b, and c are variables taking their values in the set of 
labels {1,2,3}. When a, b, and c appear together in a 
formula, they form a partition of {1, 2, 3}, so they take al- 
ways different values, and the formula is understood for 
all the different values of these variables automatically. 
(However, sometimes a formula is symmetric under the 
interchange of two such variables in which case we keep 
only one of the identical formulas.) 

The Hilbert space of a three-qubit system is H = T-L 1 <S> 
H 2 (gi'H 3 , where, after the choice of an orthonormal basis 
{|0),|1)} c H a , U a = C 2 . The \ip) e n state vectors 
are not required to be normalized in this section, so the 
€ W zero vector is also allowed. (The physical states 
arise, however, from normalized vectors.) 



A. SLOCC classification by LU-invariants 

It is a well-known and celebrated result of Diir, Vidal, 
and Cirac [9] that "three qubits can be entangled in two 
inequivalent ways." More fully, there are 1 + 1 + 3 + 1 + 
1 three-qubit SLOCC classes, that is, subsets invariant 
under LGL transformations. 

• V Nu11 (class Null): The zero- vector of H. 

• V 1 ' 2 ' 3 (class 1 1 2 1 3) : These vectors are fully separa- 
ble, which are of the form \ipi) <£> \tp2) ® \fp3)- 

« ya\bc ^hree biseparable classes a\bc), for example, 
|V>i) ® Ife) € V 1 ! 23 , where \1p23) ^ |Va) ® IV's)- For 
such | ip2s), a representative element is the standard 
B (Bell) state, 

|B) = i=(|00) + |ll». (la) 

• V w (Class W): This is the first class of genuine tri- 
partite entanglement, when no subsystem can be 
separated from the others. A representative ele- 
ment is the standard W state, 

|W) = -^(|100} + |010) + |001». (lb) 

• V GHZ (Class GHZ): This is the second class 
of genuine tripartite entanglement, the class 
of Greenberger-Horne-Zeilinger-type entanglement. 
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= 
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> 
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> 


> 



TABLE I. SLOCC classes of three-qubit state vectors identi- 
fied by the vanishing of LU-invariants Q . 

A representative element is the standard GHZ 
state, 

|GHZ) = -L(|000) + |111}). (lc) 

Formally speaking, these classes define disjoint, LGL- 
invariant subsets of H, and cover % entirely: % = 

yNull y yl|2|3 y yl|23 y V 2|13 y y3|12 y yW y yGHZ Ex _ 

cept V Nu11 , these classes are not closed. For the partial 
separability issues, we define V 123 = V w U V GHZ . 

For any £ 'H, it can be determined to which class 
belongs by the vanishing of the following quantities: 
the norm, 

n(V0 = HV-II 2 , (2a) 

the local entropies, 

s a (ip) = 4det[tr 6c (|^)<V|)], (2b) 

[here we use a normalized quantum-Tsallis entropy of pa- 
rameter 2 (see in Sec. IV Aj) . although every entropy does 
the job, since they vanish only for pure density matrices] 
and the three-tangle, 

r(V0=4|Det(V0l, (2c) 

which is given by Cayley's hyperdeterminant Det(-0) 
(l2l |H, [30j. All of these quantities are LU invariants, 
[which is U(2) x3 in this case,] moreover, n is invari- 
ant under the larger group U(8) and r under [U(l) x 

SL(2,C)] X3 U(l) x SL(2,C) x3 . It follows from the 
invariance properties and other observations [9j that the 
SLOCC classes of pure three-qubit states can be deter- 
mined by the vanishing of these quantities in the way 
which can be seen in Table |H 

Our aim is the characterization of the mixed states by 
the vanishing of some quantities, in a similar way that 
the conditions in Table H] for the quantities in ([2]) charac- 
terize the pure states. To obtain such a characterization 
scheme, we need, on the one hand, the generalization of 
classes defined somehow, and, on the other hand, a suit- 
able set of quantities which are vanishing for some classes 
determined somehow and nonvanishing for the others. 
These two issues are strongly related, and it will turn 
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out that we can define a set of quantities which suits well 
the classification given by Seevinck and Uffink @, but a 
"more complete" set of quantities suits well an extended 
but still relevant classification, which are elaborated in 

secnni 
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= 
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yl|23 
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y2|13 
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= 
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y3|12 




7^0,30 
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7^0 
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yW 


¥=0 


^0,30 


7^0 




7^0 
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yGHZ 


#o 


7^0,30 


7^0 


7^0 


7^0 




7^0 



TABLE II. SLOCC classes of three-qubit state vectors iden- 
tified by the vanishing of LSL-covariants ([3]). 



B. SLOCC classification by LSL-covariants 



the form in which we use them: 



In Borsten et. al. revealed a very elegant corre- 
spondence between the three-qubit Hilbert space H = 
C 2 ffi C 2 <g) C 2 and the FTS (Freudenthal Triple Sys- 
tem) 0JI(J) = CffiCffi Jffl J over the cubic Jordan 
algebra J = C ffi C ffi C. The fundamental point of 
this correspondence is that the automorphism group of 
this FTS is Aut[£Dl(C © C © Q] = SL(2,C) x3 , which 
is just the relevant LSL subgroup of GL(2,C) x3 , the 
LGL-group of SLOCC equivalence for three-qubit pure 
states. (This group-theoretical coincidence arises only in 
the three-qubit case.) It has been shown [ll[ that the 
vectors of different SLOCC classes of entanglement in 
the three-qubit Hilbert space are in one-to-one correspon- 
dence with the elements of different rank in the FTS. The 
rank of an element of an FTS is characterized by the van- 
ishing of some associated elements, which are covariant — 
maybe invariant — under the action of the automorphism 
group, resulting in conditions for the SLOCC classes in 
the Hilbert-space by the vanishing or non vanishing of 
SL(2, C) x3 tensors. Hence, this classification is man- 
ifestly invariant under SLOCC equivalence [111 ], which 
cannot be seen directly in the conventional classification, 
since the s a local entropies are scalars only under U(2) x3 . 
(However, the invariance of the vanishing of the functions 
s a follows easily from the fact that the local rank is in- 
variant under invertible transformations |9j.) 

Let the three-qubit state & % be expressed in the 
computational basis {\ijk) = \i) ffi \j) <8> \k)} as 



|0)= J2 4> ijk W). 



i,j,k—Q 



We can assign an element ip g 9Jt(C ffi C ffi C) to this 
and calculate some associated quantities needed for the 
identification of its rank. Here we list these quantities in 



- WW£|rf' in / ran 'r' fc (3a) 

[71 = s jf e kk ,^ k ^' k ', (3b) 

[72(V0P" =£fefe'£^V Y *\ (3c) 

[ l3 m kk ' =e w e jr ^ k r'^ k ', (3d) 

= - E mm .e nn ,e lv ^ n ^ mn 'i> im ' k (3e) 

Q\T/ &kk' E-rnm' ^nn' 



(3f) 



(For the basic definitions of Jordan algebras, Freudenthal 
triple systems and the operations and maps defined on 
them, see in [TlT | and in the references therein.) Here the 
summation for the pairs of indices occurring upstairs and 
downstairs are understood, and 
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-1 



is the matrix of the Sp(l) = SL(2)-invariant non- 
degenerate antisymmetric bilinear form: Since M t eM = 
edet(M), index contraction by e is invariant under 
SL(2, C) transformations. This shows that if we regard 
ijj and <fi as tensors that transform as a (2,2,2) un- 
der SL(2,C) x3 , then so do T (V>) and T(ip,ij;,ijj), while 
7i(V0> 72(0): an d 73 (VO; being symmetric, transform as 
(3,1,1), (1,3,1), and (1,1,3), respectively, and q{ip) 
transforms as (1,1,1); that is, it is scalar. [Note that 
for any 2x2 matrix M, the determinant 2det(M) = 
e u ,£ jj' M* J ' M*'- 7 '' , so 2 detb (?/>)] = qty).] 

The main result of [Tl| is that the conditions for the 
SLOCC classes can be formulated by the vanishing of 
these tensors in the way which can be seen in Table [TTJ 
In the light of the conditions by the norm and the four 
determinants, (see in Table [I]), this scheme constructed 
by seven quantities seems to be redundant. Indeed, it is 
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redundant for pure states, but it will turn out that this 
way leads to the generalization for mixed states. 

C. SLOCC classification by a new set of 
LU-invariants 

To follow this way, we need quantities which can be 
extended from pure states to mixed states by the convex 
roof construction There is no natural ordering 

on the tensors of ([3j so convex roof construction does not 
work directly for them, but we can form quantities from 
them taking values in the field of real numbers. During 
this, we lose the covariance under SL(2,C) x3 , but gain 
the invariance under the group U(2) x3 . 

Returning from the FTS language to the Hilbert space 
language, we have "state vectors" 

1 

i,j',fe=0 
1 

and "local operators" 
i 

7ity)e= £ biWelVI^^'leLin^ 1 ), 

i,i'=0 
1 

72 We = £ Mm^m'l e Lm(H 2 ), 
i 

7&ty)e= £ [73We]Vl^(^|eLin(^ 3 ), 

fc,fc'=0 

associated with \ip) € "H through ([3]). These are just com- 
putational auxiliaries, not state vectors and local opera- 
tors in the ordinary sense, because they depend nonlin- 
ear ly on the state vector \ip). [Note that e € H" 1 * ®rl a * = 
Un(H a -> H a *) ~ BiLin(-H a xTiMC), while 7a (^) e 
H a (8)-H a = Lin(-H a * -> V. a ), so 7a (^)e € Lin(-H a -> H a ).] 
Now, the vanishing conditions of the tensors (0 in Ta- 
ble |TT] can be reformulated. Clearly, ip = if and only if 
||t/i|| 2 = 0. Taking a look at T^ip) in (15a)) it turns out 
that T^(tp) can be written in the Hilbert space language 
as 

|T W) = YtyM 
with the "operator" 

Y(ip) = --f 1 (ifj)e<E)l<E)l-l<g)j 2 {ip)£<E)l-l<E)l<E)l3{ip)£- 
Using this, the vanishing condition of T ^{ip) for all 0, 

|T (^)) = V|0) Y(il>)\cf>) = V|0) 

^ Y(iP) = 

||y(VOI| 2 = for any norm, 



so we can eliminate the quantors and <fi from the condi- 
tion. Using the usual complex matrix 2-norm ||M|| 2 = 
tr(MtM), we have 

\\YM\\ 2 = 4(ll7iWH 2 + II72WII 2 + haWII 2 ). 

This formula has a remarkable structure, namely if we 
note that s a (ip) = hb{4>)\\ 2 + ||7cO)l| 2 and 7a(V0 = if 
and only if ||7a('0)l| 2 — 0. Now turn to the vanishing of 
T(ip, ?/>, ip), given in (l3"e|) . Again, this vanishes if and only 
if its norm ||T(^/>, if), ip)\\ 2 does. This can be calculated by 
the use of the form 

|T(^,V») = -7i(#®I®#> 
= -I(g)7 2 ('0)e®I|^> 
= -I(g)l0 73(V , )£|^> 

= \y{m)- 

(The quantity HTW: 4>, VOII 2 also appears in the twistor- 
geometric approach of three-qubit entanglement, it is 
proportional to wabc in I14I .1 About the scalar q, note 
that q(ip) — — 2Det(i/>) |l2j, and it vanishes if and only 
if the three-tangle (f2c| does. 

Summarizing the observations above, it is useful to de- 
fine the following set of real- valued functions on H.: 



n{tp) 


= U\\ 2 , 


(4a) 


2/ WO 


= |(ffiWO+saWO+0aWO)> 


(4b) 


«a(V0 


= 9b WO +9c(i>), 


(4c) 


ffa(V0 


= haWII 2 , 


(4d) 


m 


= 4||T(V,V,V0I| 2 , 


(4e) 


r 2 (V0 


= 4|<zW)| 2 . 


(4f) 



[The explicit forms of these functions and their relations 
to other important quantities can be found in the Ap- 
pendixes IA 11 IA 21 and IA31 The constant factors have 
been chosen so that < y(ip), s a (ip), g a {ip), t(if>), 7-2 W) < 
1 for normalized states, which is shown in Appendix I A 41 ] 
These quantities are obtained by index contraction of 
-0 ufc s and complex conjugated (ip l J k ) s by 8u>s from 
the tensors in ([3]) , which were obtained by index contrac- 
tion of ip l i k s and ip 1 3 ' k s by £«'S. From the contractions 
of free indices of the tensors in ([3]), we have UUU = 8 
for U £ U(2). From the contractions inside the tensors 
of ([3]), we have U t eU — edet(U) but for every factor 
det(U) there is a conjugated det(t/)* = 1/ det(U) from 
U^eU* = edet(U)* . Consequently, all the functions in 
(UJ are LU invariant, while their vanishings are still LSL 
invariant. (Again, n is invariant under the larger group 
U(8), and r 2 under [U(l) x SL(2,C)] x3 .) 

Now, the conditions for the SLOCC classes by the van- 
ishing of the tensors in ([3j (see in Table [TTJ) can be refor- 
mulated by the vanishing of the functions in Q in the 
way which can be seen in Table Hill We call the functions 
in Q pure state indicator functions for the three-qubit 
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TABLE III. SLOCC classes of three-qubit state vectors identified by the vanishing of the pure-state indicator functions given 
in ©. 



case. We will give the exact definition of indicator func- 
tions for the general case later (in Sec. IVI Cj) . until that 
point we just use this name for non-negative functions 
having the vanishing properties given in Table IIIII Al- 
though this scheme constructed by ten quantities is even 
more redundant than the previous two, but it will turn 
out that these ten indicator functions (UJ) will be neces- 
sary in the case of mixed states. Moreover, investigating 
Table Hill we can catch all the ideas leading to the general 
construction. 



III. MIXED THREE-QUBIT STATES 

Here we recall and extend the PSS classification for 
three qubits. The main concept here, first given in @, [1], 
then used and extended in [l(|, is that we define a 
density matrix to be the element of a class according to 
whether it can or cannot be mixed by the use of pure 
states of some given kinds. 



A. Convex subsets 

Let us introduce some convenient notations. The set 
of states T> = T>(H) C Lin("H) is the convex body of 
positive semidefinite operators of unit trace acting on 
while the set of pure states V C V is the set of extremal 
points of T>, which are the projection operators of rank I. 
(For the sake of simplicity, we have restricted ourselves to 
the operators of unit trace, that is, density matrices, in 
spite of the fact that the construction could be extended 
for the whole cone of positive semidefinite operators.) 
Disjoint subsets in V given by unit vectors of different 
SLOCC classes are 



properties are considered, define 



pl|2|3 = 


{MM 


M e v 1 ! 2 ' 3 


IH| 2 = 1}, 


(5a) 


■jya\bc 


{MM 


M e v a|bc , 


W 2 = i}, 


(5b) 


pw = 


{MM 


M g v w , | 


^ll 2 = l}, 


(5c) 


•pGHZ _ 


{MM 


IV.) G V GHZ 


IH| 2 = 1}: 


(5d) 



V 
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{MM I Me V 123 , IMI 2 = i}, (5e) 



which cover V entirely: V = •p 1 l 2 l 3 u'P 1|23 U7 ,2|13 U7' 3|12 U 
pw y -pGHZ gggidgg these, if only partial separability 



so V 123 = V w U p GHZ . Except P 1 ! 2 ! 3 , none of the above 
sets are closed. 

The notion of /c-separability and afe-separability [(|, 
and the relevant classes of [lfj] for three-qubit systems 
can be formulated as the convex hulls of some of the 
sets ([5]). The ^-separable states (T) 3 ~ scp ), or, equivalently 
I|2|3-separable states (I? 1 ! 2 ! 3 ) can be mixed from the 
pure states of "P 1 ' 2 ' 3 , i.e., they are fully separable. The 
a\bc-separable states (T> a \ bc ) can be written in the form 
EiPiftM ® Qbci, [Qa,i G V(n a ), Q bc ,i e V{H b ®H c )l 
where we demand only the split between a and 6c, but 
s split between b and c can also occur in the pure- 
state decompositions, so they can be mixed from the 
pure states of P 1 ! 2 ! 3 and p a l &c . The 2-separable states, 
also called biseparable states (2? 2_sop ) are of the form 
TliiPiQai,i ® QbiCi.ii so they can be mixed from the pure 
states of P 1 ! 2 ' 3 , P 1 ! 23 , P 2 ! 13 , and V^ 12 . These states are 
also of relevance because, although they are not separable 
under any a\bc split, there is no need of genuine three- 
qubit entangled pure state to mix them [6|. From the 
point of view of convex hulls of extremal points, it can 
be seen better than originally in 0] that we can define 
three new partial separability sets "between" the a | fee- 
separable and 2-separable ones. For example, the 2 1 13- 
3\12-separable states (2? 2 I 13 > 3 I 12 ) are the states which can 
be mixed from the pure states of P 1 ! 2 ! 3 , P 2 I 13 , and "P 3112 . 
States of this kind are also of relevance, since there is no 
need of l|23-separable pure states to mix them, that is, 
entanglement between the 2 and the 3 subsystems. Be- 
yond these, we use the set of W-states 

(£> W ) which 

can be expressed as the mixture of the pure states P 1 ! 2 ! 3 , 

pl|23 j p2|13 ; p3|12 ) and pW^ g0 there ig nQ need foj . pm . e 

states of GHZ type to mix them, and the set of GHZ 
states (X> GHZ ) or, equivalently, 1-separable (V 1 ^?), or 
123-separable states (I? 123 ), which is equal to the full set 
of states (T>). Summarizing, we have the following PSS 
subsets in T> arising as convex hulls of pure states of given 
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2,123 
2,2-scp 



j)2|13,3|12 2?1|23,3|12 pl|23,2|13 




1123 



2)2113 V 3 
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Pl|2|3 

FIG. 1. Inclusion hierarchy of the PSS sets T> " given in ((6]). 
kinds: 



2,11213 


= Conv^ 11213 


) =x> 3 - s °p, 


(6a) 




= Conv^ 11213 


u-p a|bc ), 


(6b) 


q-yb\ac,c\ab 


= Conv(7> 11213 


U 2? fc l ac u 2? c l ab ) ; 


(6c) 


•^2-sep 


= Conv^ 11213 


y p 1 ! 23 y p 2 ! 13 y 2? 3 I 12 ^ ; 


(6d) 


v w 


= Conv(7> 11213 


y p 1 ! 23 y p 2 ! 13 y 2? 3 I 12 


(6e) 



V 



123 



= Conv^ 11213 U P 1123 U P 2113 U V 3112 
y 2? w y 2? GHZ ) = 2? 1 - sop = £>. 

•pl23 



(6f) 



These sets are convex and they contain each other in a 
hierarchic way, which is illustrated in Fig. [TJ 

From an abstract point of view, we form the convex 
hulls of closed sets [10|, and the convex hulls of all the 
possible closed sets arising from the unions of the V" sets 
([5]) of extremal points are listed in above. We mean 
the PSS classification involving the PSS subsets (|5a]) - (f6fj) 
[and the PS classification involving the PS subsets (|6a[) 
(|6d[) and (|6fj) ] to be complete in this sense. As special, 
noncomplete cases, we get back the classification involv- 
ing only the sets T> k ~ sop and T> ak (for any /c-partite split 
ctk) obtained by Seevinck and Uffink [§[, the classification 
involving only the sets T> ak obtained by Diir and Cirac 
0, HI and also the classification involving only the sets 
jykscp anc j 2? w , obtained by Acin, BruB, Lewenstein and 
Sanpera [Io| . 



B. Classes 

Now, we determine the PSS classes of three-qubit 
mixed states. The abstract definition of these classes 
[f| is that they are the possible nontrivial intersections 
of the X>" convex subsets listed in ([6]). Since we want 
to deal also with the sets x> b \ ac ' c \ ab , we cannot draw an 
expressive "onionlike" figure as was done in Q for the 
sets V 1 ^ 3 , V a \ bc , and X> 2 - scp . We have to proceed in a 
formal manner. 

If we have the sets A\, A2, . . . , A n , all of their possible 
intersections can be constructed as the intersections for 
all i the set Aj or its complement A{. We have 10 PSS 
subsets 2>", so we can formally write 2 10 = 1024 possible 
intersections in this way. If B C A, then B n A = 0, so 
some intersections will be automatically empty ( "empty 
by construction") and, using the inclusion hierarchy of 
PSS subsets in Fig. [TJ we write only the intersections 
which are "not empty by construction." The number of 
these will turn out to be only 21. (Again, if B C A, then 
B n A = B and B n A = A, so we can write these 21 
classes as intersection sequences much shorter than 10 
terms.) Since the appearance of the D b l ac ' c l ab -type sets 
in the intersections makes the meaning of the classes a 
little bit involved, we write out the list of the PSS classes 
with detailed explanations. 

First, the class 



C 3 = £>1| 2 I 3 



(7a) 



is the set of fully separable states. 

Then come the 18 classes of 2-separable entangled 
states, that is, the subsets in V 2 ~ scp \ V 1 ^ 3 . The first 
one of them is 



C 2.s = glials n p 1 ! 23 n 2? 2113 n £> a|12 
= (x> 1123 n x> 2 ' 13 n x> 3112 ) \ x> 11213 , 



(7b) 



which is the set of states which can be written as In- 
separable states (i.e., convex combinations of "P 1 ' 2 ' 3 and 
"P 1 ! 23 pure states; the formation is not unique) and can 
also be written as 2|13-separable states and can also be 
written as 3|12-separable states but cannot be written as 
l|2|3-separable states. The existence of such states was 
counterintuitive because, for pure states, if a tripartite 
pure state is separable under any a\bc bipartition, then 
it is fully separable. For mixed states, however, explicit 
examples can be constructed [lol l3l| , which can be writ- 
ten in the form ^2 i PiQ a ,i (8 Qbc,i for any a\bc bipartition, 
but cannot be written in the form X^Pifi'M ® ® 03,i- 
Alternatively, we can say that states of this class can not 
be mixed without the use of bipartite entanglement, but 
they can be mixed by the use of bipartite entanglement 
inside only one bipartite subsystem; it does not matter 
which one. (This is class 2.8 in [6j.) The next three 
classes are 



q2. 7. a _ q~)a\bc r-i jyb\ac p. jyc\ab 

_ (j)b\ac p| q-)C.\ab^ ^ j^a\bc ^ 



(7c) 
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For example, C 211 is the set of states which can be 
written as 2|13-separable states and can also be written 
as 3|12-separable states but cannot be written as In- 
separable states. Alternatively, we can say that states of 
this class cannot be mixed by the use of bipartite entan- 
glement only inside the 23 subsystem, but they can be 
mixed by the use of bipartite entanglement inside either 
the 12 or the 13 subsystems; both of them are equally 
suitable. (These three classes are classes 2.7, 2.6, and 2.5 
in Q.) The next three classes are 



q2. 6. a _ jya\bc p q-)b\ac p J)c\ab pi jyb\ac,c\ab 



(7d) 



= J) a \ bc p| \pb\ac,c\ab ^ (pb\ac y £)c|ab)j ^ 

For example, C 261 is the set of states which can be 
written as l|23-separable states and can also be writ- 
ten as states of a new kind: where the state can be 
written as 2|13-3|12-separable states which are neither 
2|13-separable nor 3|12-separable. And this is the nov- 
elty here. Alternatively, we can say that to mix a state 
of this class we need bipartite entanglement either inside 
the 23 subsystem, or inside both of the 12 and the 13 sub- 
systems. (The latter seems like a roundabout connecting 
the 2 and 3 subsystems through the 1 subsystem.) The 
next three classes are 



£>2.5.a _ j~)a\bc p q~)b\ac p J)c\ab p J)b\ac,c\ab 



spa\bc p q~)b\ac,c\ab — ^a\bc \ 



b\ac,c\ab 



(7e) 



For example, C 2 ' 51 is the set of states which can be writ- 
ten as l|23-separable states but cannot be written as 
2|13-3|12-separable states. Alternatively, we can say that 
states of this class cannot be mixed by the use of bipartite 
entanglement only inside both of the 12 and 13 subsys- 
tems, contrary to C 2 61 . (The roundabout does not exist 
here.) (The unions C 26a U C 2 - 5a = V a \ hc n V h \ ac n V c \ ah 
are classes 2.4, 2.3, and 2.2 in [gj].) The next class is 



C 2.4 =X ,l|23pX>2|13pX>3|12 

p 2}2|13,3|12 p pl|23,3|12 p pl|23,2|13 
_/.£>2|13,3|12 p pl|23,3|12 p £>1|23,2|13\ 

\(^ 23 UP 2 I 13 UP 3 I 12 ) (7f) 

_ p£,2|13,3|12 ^ ^£>2|13 y £>3|12)j 
p rpl|23,3|12 ^ ^pl|23 y £>3|12\1 
p ^l|23,2|13^pl|23 UI) 2|13^ ) 

which is the set of states which can be mixed by the 
use of bipartite entanglement inside any two bipartite 
subsystems, but cannot be mixed by the use of bipartite 
entanglement inside only one bipartite subsystem. The 
next three classes are 



£>2.3.a _jy a \bc p J)b\ac,c\ab p jya\bc,c\ab p j)a\bc,b\ac 
|^2}&|frc,e|a& ^ ^q^c\ab y q^a\bc^ 

p |^a|frc,fc|ac ^ ^j^a\bc y ^)^l ac ^j ^ -j^b\ac,c\ab 



For example, C 2 31 is the set of states which can be mixed 
by the use of bipartite entanglement inside the 23 sub- 
system together with bipartite entanglement inside either 
the 12 or the 13 subsystems, but cannot be mixed by the 
use of bipartite entanglement inside the 12 and the 13 
subsystems only. (Note that mixing by the use of only 
one kind of bipartite entanglement has already been ex- 
cluded.) The next three classes are 



(7h) 



q2. 2. a _ jyb\ac,c\ab p jy a \bc,c\ab p J)a\bc,b\ac 
qyb\ac,c\ab ^ ^q^a\bc.c\ab y q^a\bc.b\ac\^ 

For example, C 231 is the set of states which can be mixed 
by the use of bipartite entanglement inside both the 12 
and the 13 subsystems together, but cannot be mixed by 
the use of bipartite entanglement inside the 23 subsystem 
together with bipartite entanglement inside only one of 
the 12 or the 13 subsystems. The next class is 



(7i) 



C 2 ' 1 = £>2|13,3|12 p £>1|23,3|12 p £>1|23,2|13 p £>2-sep 

_ 2-,2-scp ^ ^-p2|13,3|12 y -pl|23,3|12 y £,1123,2113) 

which is the set of states which can be mixed by the use 
of bipartite entanglement inside all the three bipartite 
subsystems, but cannot be mixed by the use of bipartite 
entanglement inside only two (or one) bipartite subsys- 
tems. (The union C 2A UC r 



uc 2 - 3 - 2 uc 



2.3.3 



uc 2 



u 



£2.2.2 y £2.2.3 y C 2.1 = ^2-scp \ ^1|23 y £,2|13 y £,3|12) fe 

class 2.1 in [6j.) 

Then come the 2 classes of states containing genuine 
tripartite entanglement [To[, that is, the subsets in T> \ 
V 2 - se P. The class 



(7j) 



is the set of states which cannot be mixed without the 
use of some tripartite entangled pure states, but there is 
no need for GHZ type entanglement [Ioj |. The class 



C 



GHZ 



= v w n x> 3 = v 123 \ v 



VW 



(7k) 



is the set of states which cannot be mixed without the 
use of GHZ type entanglement. All the above classes are 
PSS classes. For the PS classification define the class of 
states containing genuine tripartite entanglement instead 
ofC w andC GHZ : 



C 1 = C w U C GHZ = I? 123 \ 2? 2 - sc P 



(71) 



(7g) 



Except C 3 , the C " PS(S) classes above are neither 
convex nor closed, but, by construction, they cover V 
entirely. Unfortunately, we cannot draw an onionlike fig- 
ure illustrating these classes, like the one in 6j (maybe it 
could be drawn in three dimensions); we only summarize 
these 1 + 18 + 1 + 1 classes in Table HVl The non-emptiness 
of the PS(S) classes above is not obvious, since it depends 
on the arrangement of different kinds of extremal points. 
(We know only that they are not empty by construction.) 
This issue has not been handled yet, but experiences in 
the geometry of mixed states [l| suggest that the arrange- 
ment of different kinds of extremal points leading to some 
empty classes would be very implausible. 
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TABLE IV. PSS classes of mixed three-qubit states and PS classes of mixed tripartite states. Additionally, we show the 
classifications obtained by Seevinck and Uffink [f|, Diir, Cirac and Tarrach and Acfn, Brufi, Lewenstein and Sanpera [Io| . 



C. Convex roof quantities 

As a next step, we obtain indicator functions on mixed 
states from the pure-state indicator functions (j4|) by con- 
vex roof construction p~5l - fT7| . In general, let 



/ :V — > M 

be a continuous function. Then its convex roof extension 
is defined as 



V 



min 2^ pi f(i>i), 



(8) 



where the minimization takes place on all pure-state de- 
compositions of g: < pi, J2iPi = 1) X^I^X^I = 8- 
The existence of the minimum is crucial for our construc- 
tion. It follows from the Schrodinger mixture theorem 
[32| . also known as the Gisin-Hughston-Jozsa-Wootters 
lemma [U [34| , that the decompositions for m pure states 
are labeled by the elements of the compact complex mani- 
fold, called Stiefel manifold, St d (C m ) = U(m)/U(m-d), 
where d = dim'H The Caratheodory theorem en- 
sures that we need only finite m. or to be more precise 
m < d 2 , shown by Uhlmann [35|]. These observations 
guarantee the existence of the minimum in (J8J . 

Now, it is easy to prove the following necessary and 
sufficient conditions for the PSS subsets ^ given by the 
convex roof extension of the indicator functions (U): 





y u (e) 


= o, 


(9a) 






= o, 


(9b) 


g £ j)b\ac,c\ab ^_ 


=> 9a (e) 


= 0, 


(9c) 


g G £> 2 ~ sop <S= 


=> t u (e) 


= 0, 


(9d) 


gev w <= 




= 0. 


(9e) 



To see the =>■ implications, observe that all the 2> ' ■ PSS 
subsets are the convex hulls of such pure states [see in 



©] for which the given functions vanish [see in Table 
IIII| . Since these functions can take only non- negative val- 
ues, the minimum in the convex roof extension is zero. 
To see the <= implications, note that if the convex roof 
extension of a non-negative function vanishes then there 
exists a decomposition for pure states for which the func- 
tion vanishes. Again, the vanishing of a given function 
singles out the pure states [see in Table HIT] from which 
the states of the given 2> ■ ■ PSS subset can be mixed [see 
in ©]. 

The necessary and sufficient conditions for the PSS 
subsets © yields necessary and sufficient conditions for 
the PSS classes, and we can fill out Table fVl for the iden- 
tification of the PSS classes of Table UVT given for mixed 
states, similar to Table iHfl given for pure states. Because 
of their vanishing properties, we call the convex roof ex- 
tension of pure indicator functions mixed indicator func- 
tions. 

Note that the convex roof extension is a nonlinear op- 
eration: (fi + / 2 ) u 7^ /i u + /2 U - However, an inequality 
holds, for example, s a = g\, + g c and s^ 1 = (<?& + g c ) u > 



■g^, so can be nonzero even if both g^ and g^ are 



zero. This is why we could identify 21 classes of mixed 
states by the use of the convex roof extension of func- 
tions which identify only 6 classes of state vectors. On 
the other hand, if a classification does not involve all the 
PS(S) subsets, then, through ©, we have to use only 
some of the indicator functions, for example, y, s a and t 
for the classification obtained by Seevinck and Uffink [6] , 
y and s a for the classification obtained by Diir, Cirac and 
Tarrach and y, t and t 2 for the classification obtained 
by Acm, Brufi, Lewenstein and Sanpera fioj j . 



IV. EXAMPLES 

At this point, the most important question is whether 
all of the PSS classes in (J7J are nonempty. Of course, this 
can be checked by the use of ©, but calculating convex 
roof extensions symbolically is a hard problem. Here we 
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= 
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> 


> 


> 


> 



TABLE V. PSS classes of mixed three-qubit states given in table ITVl identified by the vanishing of the mixed indicator functions 
(convex roof extension of the indicator functions ((Jl). 



give considerations apart from convex roofs. 

The classes given by Seevinck and Uffink in [|| are 
nonempty, which are C 3 , C 2 8 , C 2 ' 7 a , the unions C 26a U 
C 2ba , the union C 2A U C 231 U C 232 U C 233 U C 221 U 
C 2. 2 .2 uC 2.2.3 uC 2.i j an dC w UC GHZ . Both of the classes 

C w and C GHZ are nonempty [l(|. On the other hand, 
the pure sets ([SJ are contained by the following classes: 

pl\2\3 c £3 jya\bc ^ q2. 5. a -pW qW -pGHZ <- £>GHZ 

so we have additionally that C 2 5 a is nonempty. In 
the next paragraphs, we construct states contained in 
classes C 2 2 a and C 21 . This justifies the use of b\ac-c\ab- 
separable sets in the classification (since we can distin- 
guish between C 2 2 a and C 21 by the use of these), al- 
though the nonemptiness of C 2S a , C 2A , and C 23a has 
not been shown yet. 

From the point of view of "mixtures of extremal 
points," it is easy to check that the bipartite subsystems 
are separable for states in some PS subsets as follows: 



g G £> a|bc 

g qjb\ac,c\ab 



f?23 sep. and gi3 sep. and Q12 sep. 

Q ac sep. and g ab sep. 

Qbc sep. 



Unfortunately, the reverse implications are not true. For 
example, for the standard GHZ state ((Tell , all bipar- 
tite subsystems are separable, although |GHZ)(GHZ| £ 
I? 1 ! 2 ! 3 . However, the negation of the implications above 
will turn out to be useful: 



Q £V^ 3 

p 4 J) h \ ac A ah 



£>23 ent. or ^13 ent. or gy^ ent. 

g ac ent. or g ab ent. 

Qbc ent. 



The entanglement of two-qubit subsystems can be easily 
checked, for example, by the Peres-Horodecki criterion 
3G, 371: 



lo separable 



> 0. 



(10) 



Here w G T>(H b ®H C ), and tj means transposition on the 
first subsystem, which, although it is basis-dependent, 



the positivity of a;' 1 is not. The <= implication in (TTTJ)) 
holds only for qubit-qubit or qubit-qutrit systems. 

Now, take a g G 2? 2 I 13 ' 3 I 12 . Then g 2 3 is always sep- 
arable, but if both Q12, and 013 are entangled, then by 
the above observations we have g 4_ I? 1 ! 23 , g 4_ I? 2 ! 13 , 
q £ P 3 ! 12 , moreover, g 4_ V 1 ^ 13 , and g 4_ V x \ 23 ^ 12 . 
This singles out exactly one class from Table ITVl namely 
C 2 - 21 . So if we can mix a state g from V 1 ^ 3 , V b \ ac , 
and •p c l ab ; whose g a b and g ac subsystems are entangled, 
then g G C 22a . For example, such a state is the uniform 
mixture of projectors to the |0)f, <g> |B) ac and |0) c (g> |B) j 
vectors: 



i|0)(0| c ®|B)(B| ah eC 2 - 2 - a , 



||0><0| 6 03 1B)<B| 



where |B) is the usual Bell state ((la)) . 

Now, take a g G V 2 ~ scp . Then if the states of all the 
two-qubit subsystems are entangled, by the above obser- 
vations we have g 4_ I? 1 ! 23 , g 4_ I? 2 ! 13 , g 4_ I? 3 ! 12 , more- 
over, q i V 2 ^ 3 ^ 2 , g £ V 1 ^ 12 , and g £ V 1 ^ 13 . 
This singles out exactly one class from Table ITVl namely 
C 21 . So if we can mix a state g from P 1 ! 2 ! 3 , "P 1 ! 23 , "P 2 ! 13 , 
and P 3 ' 12 , whose all two-qubit subsystems are entangled, 
then g G C 21 . For example, such a state is the mixture 
of projectors to the previous two vectors together with 
|1)„® |B) 6c : 



i|0)<0| 6 ®|B)<B| ac + j 



:®\B)(B\ ab 



-IIXIU^IBXB^GC 2 ' 1 . 



GENERALIZATIONS I. 
SUBSYSTEMS 



THREE 



The considerations written out in detail in Sees. iHl and 
IIIII contain the main ideas which will be generalized in 
this and in the next sections. In this section we break 
up with qubits, and consider tripartite systems composed 
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from subsystems of arbitrary dimensions. Obviously, this 
has no influence on the PS sets and PS classes, given in 
Sees. IIII Al and UlI B| the only question is about the con- 
struction of mixed-state indicator functions of Sec. IIII CI 
The generalization to arbitrary number of subsystems is 
left to the next section. 



A. Pure state indicator functions for tripartite 
systems from the FTS approach 

To get the necessary and sufficient conditions for the 
PS classes in the tripartite case, we need the generaliza- 
tions of the pure-state indicator functions in (|4bp - (l4e|) . 
Apart from continuity, the main — and only — requirement 
for these is to satisfy the vanishing requirements for 
pure states given in Table IIIII (apart from the column 
t 2 , and for row V 123 instead of rows V w and V GHZ ). 
Then their convex roof extensions satisfy the vanishing 
requirements for mixed states given in Table |V] (apart 
from the column t 2 , and for row C 1 instead of rows 
C w and C GHZ ), since in ([5a ]) -([9d | we have used only the 
vanishing-requirements for pure states. 

The pure-state indicator functions of §4§ have been ob- 
tained in the FTS approach, which works only for the 
qubit case. However, some parts of the definitions can be 
generalized. To do this, our basic quantities will be the 
local entropies s a (ip) — S q (jt a ) instead of the functions 
g a (ifi) given in (|4dl) , since the former ones are defined for 
all dimensions. [Here we use the notation for the pro- 
jector 7r = and 7r a = tr;, c (7r).] The most basic 
quantum entropy is the von Neumann entropy, 



S(g) = -tr [filing)], 



(11) 



having the strongest properties among all entropies. The 
Tsallis entropy, sometimes called g-entropy, in the quan- 
tum case is defined as 



S, 



1 



q ( e ) = — [tr( e «)-l], g>0, (12a) 



which is a nonadditive generalization of the von Neumann 
entropy: lim q -yi S q (g) = S(g). Again, as in (|2bj) . we 
can use the concurrence-squared, which is the normalized 
Tsallis entropy of parameter 2: 



= ^ri^(e) = j-j [i - 1 V)] , (12b) 



if we prefer to deal with polynomials in the tp^ and 
(i/j l 3 k ) coefficients. This is the nontrivial polynomial 
of the lowest degree which is also an entropy, that is, 
Schur-concave, so tells us something about mixedness. 
[In (|12bj) . d is the dimension of the Hilbert space on which 
g acts, so the prefactor -j4t ensures that < C 2 (g) < 1.] 
Obviously, for all Tsallis entropies of the subsystems, 
s a(V0 = S q (Tv a ) fulfils the corresponding column of Table 
IIIII since it vanishes if and only if the subsystem is pure, 
which means the separability of that subsystem from the 



rest of the system if the whole system is in pure state. 
From Pc|) and (|4d)) . it turns out that y, given in (|4b|) . is 
just the average of the local entropies y = g(si + s 2 + 
S3), vanishing if and only if no entanglement is present. 
This works well not only for qubits, so we can keep this 
definition of y. 

The functions g a in (|4dj) can also be expressed by the 
local entropies (|4"c| for qubits as g a — h(sb + s c — s a ). 
Can this definition be kept for subsystems of arbitrary 
dimensions? For V 1 ' 2 ' 3 , obviously g a = 0. For y a ' fcc , 
the subsystem a can be separated from the others so the 
subsystems a and be are in pure states, s a — and s& = 
s c =/= 0, from which g a ^ and gb = g c = 0. So the first 
five rows of the g a columns of Table IIIII is fulfilled. For 
the last row, we need that g a > when genuine tripartite 
entanglement is present. This is the problematic point. 
This question can be traced back to the subadditivity of 
the Tsallis entropies. Raggio's conjecture [38[ about that 
is twofold: For q > 1, 



S q (g) < Sqigx) + S q (g 2 ), 
S q {g) = S q {g!) + S q (g 2 ) 



Q = Ql ® Q2, 

gi or g 2 pure. 



(13a) 
(13b) 



[Note that for < q < 1, there is no definite relation 
between S q (g) and S q (gi) + S (g 2 ).] Both statements 
hold for the classical scenario [38|], which can be mod- 
eled in the quantum scenario by density matrices being 
LU equivalent to diagonal ones. The first part (| 13al) of 
the conjecture has been proven by Audenaert [39(. This 
guarantees the non-negativity of the functions g a : For 
pure states, S q (r a ) = S q (ix bc ) < S q (ir b ) + S q (TT c ), so 
< !(s& + s c — s a ) = g a . On the other hand, (|13bl) is 
exactly what we need: \ip) £ V 123 if and only if neither of 
its subsystems are pure, which means that there is subad- 
ditivity in a strict sense, so < ^(sf, + s c — s a ) — g a . The 
<= implication in (| 13b|) holds obviously, but the whole 
second part (|13b|) of the conjecture, to our knowledge, 
has not been proven yet. A very little side result of our 
work is that Raggio's conjecture holds for the very re- 
stricted case of two-qubit mixed states which are, at the 
most, of rank 2. 

We note that the von Neumann entropies (q — > 1) of 
the subsystems are not suitable for the role of the func- 
tions s a , if we want to write the functions g a by that as 
5 (sf,+s c — s a ), since the von Neumann entropy is additive 
for product states without any reference to the purity of 
the subsystems: 



S(g)<S( gi ) + S(g 2 ), 
S(g) = S( Q1 ) + S(g 2 ) 



g = gi <8> g 2 . 



(14a) 
(14b) 



Indeed, it is easy to construct a tripartite state, which is 
not separable under any partition, but has vanishing g a 
(defined by the von Neumann entropy). For example, let 
dirnH" = 4, then for the state 

M = i(|000) + |101) + |210) + |311>) 
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7T23 = 7T2 ® 7T3, SO ffl(^) = 5 (£(^2) + S(7T 3 ) - S^l)) = 

0, while S(tti) = In 4, and S(ir 2 ) = S(w 3 ) = In 2, so 
neither of the subsystems are pure, the state is genuinely 
tripartite entangled. 

The Renyi entropy is defined as 

Sf( Q ) = Y ^]n[tT( e «)}, q>0, (15) 

which is another generalization of the von Neumann en- 
tropy: lim 9 _>i Sq (g) = S(g), having the advantage of 
additivity: 

Sf(g) = Sf( gi ) + Sf(Q 2 ) Q = ei ®Q 2 . (16) 

This is an advantage when entanglement is studied in the 
asymptotic regime, when the state is present in multiple 
copies and properties are investigated against the number 
of copies. Again, this advantage is a disadvantage from 
our point of view, the Renyi entropies of the subsystems 
are not suitable for the role of the functions s a , if we 
want to write the functions g a by that as \{sb + s c — 
s a ). Moreover, subadditivity does not hold for Renyi 
entropy, so the non-negativity of the functions g a defined 
by Renyi entropies is not even guaranteed. (For further 
properties and references on the quantum entropies, see, 
for example, QUO-SI-) 



B. Pure state indicator functions for tripartite 
systems outside the FTS approach 

Fortunately, it is easy to define the pure-state indica- 
tor functions of three subsystems of arbitrary dimensions 
without the issues of equality in the subadditivity of q- 
entropies. Again, the basic quantities are the local en- 
tropies, and we use a "multiplicative" definition for the 
functions g a instead of the "additive" one, which came 
from the FTS approach, 

y(tl>) = sity) + s 2 (il>) + s 3 (fl>), (17a) 

s o (V0 = S q (ir a ), (17b) 

ffo(V0 = s b (ip)sc{ip), (17c) 

t(V) = aiWO^MVO- (17d) 

These functions obviously reproduce the relevant part of 
Table UTTl so, by (|9a|) — (|9d|) . their convex roof extensions 
reproduce Table IVl for the identification of the PS classes 
of the tripartite case given in Table IIV1 The structure 
of the formulas above give us a hint for the generaliza- 
tion for arbitrary number of subsystems of arbitrary di- 
mensions: Te just have to play a game with statements 
like "being zero," with the logical connectives "and" and 
"or," parallel to the addition and multiplication, and also 
parallel to the set-theoretical inclusion, union, and inter- 
section. 



VI. GENERALIZATIONS II. - PARTIAL 
SEPARABILITY OF MULTIPARTITE SYSTEMS 

In the previous sections, we have followed a didactic 
treatment in order to illustrate the main concept; now 
it is high time to turn to abstract definitions to handle 
the PS classification and criteria for arbitrary number of 
subsystems of arbitrary dimensions. 

For n subsystems, the set of the labels of the subsys- 
tems is L = {1,2, . . . , n}. Let a = L\\L 2 \ ■ ■ ■ |£& denote 
a /c-partite split, that is, a partition of the labels into k 
disjoint non-empty sets L r , where L\ U L 2 U ■ ■ ■ U Lk = L. 
For two partitions, ft and a, ft is contained in a, denoted 
as ft -< a, if a can be obtained from ft by joining some — 
maybe neither — of the parts of ft. This defines a partial 
order on the partitions. [It is easy to see from the defi- 
nition that a -< a (rcflcxivity) ; if 7 -< ft and ft -< a then 
7 X a (transitivity); if ft ^ a and a ^ /3 then a = (3 
(antisymmetry).] For example, for the tripartite case 
1 1 2 1 3 -< a\bc if! 123. Since there is a greatest and a small- 
est element (the full n-partite split and the trivial parti- 
tion without split, respectively, 1|2| . . . \n ^ a X 12 . . . u,) 
the set of partitions of L for ^ forms a bounded lattice. 

A. PS subsets in general 

The first point is the generalization of the PS subsets 
T)- - . Let V a be the set of pure states which are separable 
under the partition a — L\\L 2 \ . ■ . |ifc, but not separable 
under any ft -< a. Then the PS subset of a-separable 
states is 

V a = Conv (J T" 3 , (18a) 

which is a special case of the PS subsets of a- separable 
states 

V a = Conv (J |J V 13 = Conv (J V a , (18b) 

with the label a being an arbitrary set of partitions. [In 
the writing we omit the {. . . } set brackets, as was seen in, 
e.g. . (joe)) ,] The set of fc-separable states T> k ~ sep arises as a 
special case where the a elements of a are all the possible 
fc-partite splits. Note that in general, the a partitions are 
not required to be fc-partite splits for the same k. This 
freedom can not be seen in the case of three subsystems. 

The V a sets are not closed if and only if a is not the 
full ri-partite split l|2|...|n, but Up^ a V^ is closed, so 
the sets V a are closed, and convex by construction. Note 
that different a labels can give rise to the same T> a sets; 
in other words, the a. H> T> a "labeling map" defined 
by (|18b[) is surjective but not injective. For the full PS 
classification we need all the possible different T> a sets. 
Because of the nontrivial structure of the lattice of par- 
titions, obtaining all the different PS sets is also a non- 
trivial task. We cannot provide a closed formula for that, 
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but only an algorithm. Before we do this, we need some 
constructions. 

First, observe that if j3 < a then V 13 C V a [from def- 
inition (|18a|) . and the transitivity of ;<], from which it 
follows that for the labels (3 and a, if for every f3 £ (3 
there is an a € a for which (3 < a then T>@ C T> a . 
[From definition ([T8b")l . We will prove the reverse too.] 
These observations motivate the extension of ^ from the 
partitions to the labels: 



(3 d a 



dcf. 



V/3 e /3,3a € a : (3 < a. (19) 



Note that, at this point, the relation < on the labels is 
not a partial order, only the reflexivity and the transitiv- 
ity properties hold for that. The antisymmetry property 
fails, which is the consequence of that the definition (fT9|) 
was motivated by the inclusion of the PS sets, and dif- 
ferent as can lead to the same PS set. Independently of 
this problem, which will be handled later, the following 
is true: 



(3 ■< a 



V 13 C V a . 



(20) 



For the proof, see Appendix IB 11 Again, note that the 
relations -< and C are defined on nonisomorphic sets, 
so (j2"0)l does not contradict the fact that the latter is a 
partial order while the former is not. 

The next step is to define those labels for which ■< will 
be a partial order. A label a is called a proper label if 



a ^ a'. (21) 



Va, a' € a, a^a 



On the set of proper labels, the relation -< defined in (jT9|) 
is a partial order. For the proof, see Appendix IB 21 A 
corollary is that the set of proper labels for -< forms a 
bounded lattice, its greatest and smallest elements are 
the one-element labels of full n-partite split and the triv- 
ial partition without split, respectively: l|2|...|n^a^ 
12... n. 

Is it true that every PS subset can be labeled by proper 
label? Do different proper labels lead to different PS 
subsets? In other words, is the a M> V a "labelling map" 
from the set of proper labels to the set of PS subsets an 
isomorphism? The injectivity is the ^ implication from 
the observation, that for a, (3, proper labels 



(3 = a 



V f3 = V ol 



(22) 



For the proof, see Appendix IB 31 If (3 is a label, then we 
can obtain a unique proper label from that, if we drop 
every (3 g (3 for which there is a (3' E (3 for which (3 < f3' . 
The remaining partitions form a proper label which we 
denote a, and the partitions which have been dropped 
out form a label which we denote 7. Then (3 = 0:7, 
which means the union of labels a and 7. (We omit the 
union sign too.) Our next observation is useful for this 
case. For the general labels a and 7, 



which means that when we obtain a proper label a from 
a general label f3, as was done above, both of these lead 
to the same PS subset. For the proof, see Appendix IB 41 
Since all PS subsets arise from general labels, the above 
shows that they arise also from proper labels, which is 
the surjectivity of the labeling by proper labels. 

Now we have that the set of proper labels is isomorphic 
to the set of PS subsets. The former one is much easier 
to handle. Moreover, (|2"0|) states now that the lattice of 
a proper labels with respect to the partial order ^ is 
isomorphic to the lattice of T> a PS subsets with respect 
to the partial order C. (This lattice is the generalization 
of the "inclusion hierarchy" in Fig. [TJ) To get all the PS 
subsets, we have to obtain all the proper labels. A brute- 
force method for this is to form all the (3 labels (all the 
subsets of the set of all partitions), then obtain the proper 
labels a as before ([3 = otj) and keep the different proper- 
labels obtained in this way. A much more sophisticated 
algorithm is given in Appendix [ 



B. PS classes in general 

The second point is the generalization of the PS classes 
C "' , which are the possible non-trivial intersections of 
the PS subsets 2> ' • . Constructing these requires direct 
calculations for a given n, as was done in Sec. IHI Bl 

Let us divide the set of proper labels into two disjoint 
subsets, a and f3; then all the possible intersections of 
PS subsets can be labeled by such a pair, which is called 
class label, as 



(24) 



Pep 



It can happen that C— S- = by construction, under which 
we mean that its emptiness follows from the inclusion hi- 
erarchy of PS subsets. For example, if T>P C V a for some 
a. G a and (3 € f3, then the intersection above is identi- 
cally empty. The PS(S) classes for three subsystems in 
Sec. IIIIBI were obtained by the use of this observation. 
In this general framework, this observation is formulated 
as follows: 



C 



9L,P 



Pi v 13 c y v a 

PGP a£a 

3a e a, 3f3 e (3 : V p C V c 
3a e a, 3/3 e (3 : (3 < a. 



7 -< a 



(23) 



[Equivalence (i) comes from De Morgan's law A n B = 
AUB. Equivalence (ii) comes from the identity B C 
A BC\A = B\ A = ®. Implication (Hi) comes from 
B C A => B n B' C A U A'. Equivalence (iv) is ([2"0|) .] 
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Implication (iii) is the point which makes it possible 
to formulate the emptiness of PS classes by the use of 
labels only. That is still a question whether implication 
(iii) can be replaced with a stronger one, which leads to 
a condition involving only labels again. (The problem is 
that we have no interpretations of fl and U in the lan- 
guage of labels.) Our first conjecture is that implication 
(iii) above is the strongest one which leads to a condition 
involving only labels. 

Summarizing, we have 



C 



3aea,3/3e/3 : [3 < a. (25a) 



If the right-hand side holds, then we say, according to 
the conjecture above, that C— is empty by construction. 
Since this implication is only one way, it could happen 
that C— = for such class label a, (3 for which the right- 
hand side does not hold. However, we think that this 
cannot happen: Our second conjecture is that there is an 
equivalence in H25a\) : that is, all the PS classes which are 
not empty by construction are nonempty. [This implies 
the first conjecture above, but it can still happen that 
implication (iii) can be replaced by a stronger condition, 
so the first conjecture is false. Then the (I25ap definition 
of the emptiness- by-construction changes, and the second 
conjecture concerns this new definition.] The motivation 
of this is the same as in the tripartite case, (see at the 
end of Sec. IIII Bl) . where the PS classes conjectured to be 
non-empty were obtained under the same assumptions. 

An advantage of the formulation by the labeling con- 
structions is — roughly speaking — that by the use of that 
"we have separated the algebraic and the geometric part" 
of the problem of nonemptiness of the classes. At this 
point, it seems that we have tackled all the algebraic is- 
sues of the question, and these conjectures cannot be 
proven without the investigation of the geometry of T> 1 
more precisely, the geometry of the different kinds of V a 
extremal points. 

The negation of (|25aj) leads to 



Va e a,V/3 £ §_ : f3 ^ a, (25b) 



so if we obtain all a, (3 class-labels for which 
the right-hand side of this holds non-emptiness-by- 
construction"} then we will have all the nonempty 
classes, together with some empty classes if the second 
conjecture does not hold. Because of the nontrivial struc- 
ture of the lattice of proper labels, obtaining all the class- 
labels leading to nonempty-by-construction classes is also 
a nontrivial task. The number of all the partitions of n 
grows rapidly [43|, H3] , which is only the number of the 
PS subsets of a-separability V a . So the number of all the 
PS subsets T> a grows more rapidly, and the number of all 
the PS classes C—'S. grows even more rapidly. However, 
at least, it is finite. 



C. Indicator functions in general 

The third point is the generalization of the indicator 
functions. Let F : T>(H K ) — > 1 be a continuous func- 
tion for all K C L, that is, for all — also composite — 
subsystems. The only condition on F is 

F(ff) > 0, with equality if and only if g is pure, (26) 

for example, the von Neumann entropy or any Tsallis or 
Renyi entropies are suitable. (Note that the additional 
requirements of the features of LU invariance, convexity, 
Schur-concavity, additivity, being homogeneous polyno- 
mial, etc., are only optional; they will not have been 
needed for the construction.) For all K C L subsystems, 
let the following functions on pure states be defined: 



Ik ■ V — > R) 
f K (n) = F(n K ), 



(27) 



where again, tt = and ttk — tr-^(7r), with K — 

L\K. Then, for the fc-partite split a = Li\L 2 \ ■ ■ ■ \L k , f Lr 
identifies the bipartite split L r \L r , (where L r = L \ L r ,) 
as 



(28) 



which is the consequence of (|26|) . Note that a is the 
greatest element which is smaller than L r \L r for all r. 
Then the function 



(29) 



has the ability to identify the fc-partite split a as 

/ a (7T)=0 7T6 |J V?. (30) 



All non-negative f a functions satisfying (|3Tfl) are called 
a-indicator functions, not only the ones defined in (|29p . 
The generalization of (f29| for more-than-one partitions, 
that is, for all labels, is defined as 



/«w = n 



(31) 



being the generalization of ([T7|) . It vanishes exactly for 
the convenient "P Q s 



/aW=0 



ne (J U^- 



(32) 



All non-negative f a functions satisfying (|32|) are called 
a-indicator functions, not only the ones defined in (|31l) . 
For example, the functions in Q were not constructed 
by (|3ip. but still satisfy (|3"2")l . which is equivalent to the 
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relavant part of Table IIIII for the three-qubit case. Now, 
the vanishing of their convex roof extension 



identifies the PS sets: 



(33) 



being the generalization of ©. Indeed, fa(o) = if and 
only if there exists a decomposition g = ^i-Pi 71 "* such 
that foi(^i) — for all i (f a is non- negative) , which 
means that TVi € [j aea {Jp^ a (E2)i which means that 
>>■ P". 



D. Entanglement-monotone indicator functions in 
general 

There is a possibility to choose indicator functions so 
that they obey some axioms of entanglement measures 
[45| . The most fundamental one of them is the mono- 
tonicity under LOCC [H H||. A /z : 2? ^ 1R is (nonin- 
creasing) monotone under LOCC if 



A»(A(g)) < h(q) 



(34a) 



for any LOCC transformation A, which expresses that 
entanglement can not increase by the use of local oper- 
ations and classical communication. A \i : T> — > R is 
nonincreasing on average under LOCC if 



^Pj^Qj) ^ vie). 



(34b) 



where the LOCC is constituted as A = X^'Aj's wnere 
the AjS are the parts of the LOCC realizing the out- 
comes of selective measurements, and Qj — —Aj(g) with 
Pj = tr[Aj(p)]. This latter condition is stronger than the 
former one if the function is convex: 



E 



P3Q3 



< 



(34c) 



for all ensemble {(pj,gj)}, which expresses that entan- 
glement can not increase by mixing. A fi : T> — > K 
is entanglement-monotone if (|34b[) and (I34c[) hold for 
that There is common agreement Q that LOCC- 

monotonity (|34ap is the only necessary postulate for a 
function to be an entanglement measure. However, the 
stronger condition (|34b|l is often satisfied too, and it is 
often easier to prove. This holds also for our case. 

If /1 is defined only for pure states fi : V — > R, then 
only (|34bp makes sense, whose restriction is 



(35) 



Here {(pi, fti)} is the pure ensemble generated by all the 
Kraus operators of all AjS from the input state ir. (Not 



all Hi members of the ensemble are accessible physically, 
only the outcomes of the LOCC, which are formed by 
partial mixtures of this ensemble [45j]. Mathematically, 
however, we can use the pure ensemble, which make the 
construction much more simple.) If we have such a func- 
tion a : V — >• K, (|34bp holds for its convex roof extension 
SE Si: 



Y^PiV(k 1 ) < /x(tt) 



5>M U (ft)<M U (e)- (36) 



Since the convex roof extension of a function is convex 
[Eq. l34c| (moreover, it is the largest convex function tak- 
ing the same values for pure states as the original function 
does, 35]), fJ. u (g) is also entanglement-monotone. 

Now, we construct indicator functions which are 
entanglement-monotone. (These are denoted m in con- 
trast with the general fs.) This is carried out in four 
steps. 

(i) It has also been shown in j4|| |46[ that if F : 
r D(H K ) — > M. is unitary invariant and concave, then the 
fx functions defined in (|27p are non increasing on aver- 
age for pure states, that is, obey ([3"5"j) . So let 



m K (ir) = M(it k ) 



(37) 



with M : D{T-L K ) — > R vanishing if and only if the state 
is pure, as before, but now we demand also unitary in- 
variance and concavity. The von Neumann entropy (jlip . 
the Tsallis entropies (|12ap for all q > 0, and the Renyi 
entropies (fTS"]) for all < q < 1 are known to be concave 
1], and all of them are unitary invariant. 

(ii) Clearly, the functions obeying ([331) form a cone; 
that is, their sums and multiples by non-negative real 
numbers also obey (|35p . so we can conclude that the sums 
of the functions tuk also obey (|3"5|) . Here, instead of the 
original sums in (|29l) . we introduce the arithmetic mean 
of the mi r functions, 



1 fe 

m a {Tr) = -j^2m Lr (K), 

r=l 



(38) 



which are also indicator functions, since they obey (|30p . 
[The factor 1/k is not really important, but the next step, 
and in the three-qubit case y = l/3(si + s 2 + s 3 ) from 
(|4bp - (|4"c|) motivate the use of mean values.] 

(iii) The only problem we face here is that the set of 
functions obeying (|35p is not closed under multiplication, 
which is the case of the f a functions of ([31]) . This is re- 
lated to the fact that the product of two concave func- 
tions is not concave in general. Moreover, a recent result 
of Eltschka et. al. suggests that functions obeying (|35|) 
cannot be of arbitrary high degree (see Theorem I in 
471 ] . concerning a special class of functions), so we make 
a trial of such a combination which does not change the 
degree, but still fulfils the conditions (|32|) . The geometric 
mean will be proven to be suitable, which is just a root 
of the product given in (|3"Tj) 



[n m "( 7r ) 



(39) 
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where I = \a\, the number of m a s in the product. These 
functions obviously obey ([3"2j) and also (f3"5j) . which latter 
is proven in Appendix IB 61 

(iv) Now, the function m a (Tr) of (|39[) is nonincreasing 
on average for pure states (|35l) so 

m^io) = min^p. i TO a (7r 4 ) 

i 

is also non increasing on average (|34bl) [because of ([36]) ] , 
so it is entanglement-monotone and also identifies the PS 
subsets 

m^(g) = geV a , (40) 

as in (j3"3")) . 

VII. SUMMARY AND REMARKS 

In this paper we have constructed the complete PS 
classification of multipartite quantum states by the PS 
classes arising from the PS subsets (|18bl) . together with 
necessary and sufficient conditions for the identification 
of the PS classes through the necessary and sufficient 
conditions for the identification of the PS subsets (|3"3")) 
by indicator functions arising as convex roof extensions 
of the pure-state indicator functions (1511) . The indicator 
functions can be constructed so as to be entanglement- 
monotone (Sec. I VI Dp . A side track is the PSS classifi- 
cation of three-qubit states, fSec. IIIip . where a different 
set of indicator functions has been obtained CP, © by 
the use of the FTS approach of three-qubit pure-state 
entanglement. 

Now, we list some remarks and open questions, first 
about the general case, (i) As was mentioned before, 
this PS classification scheme is an extension of the classi- 
fication based on fc-separability and afc-separability given 
by Seevinck and Uffink Q , which is the extension of the 
classification dealing only with afc-separability given by 
Diir and Cirac 0, H|. (ii) The nonemptiness of the new 
classes was only conjectured. More fully, we could not 
give necessary and sufficient condition for the nonempti- 
ness of the PS classes in the purely algebraic language 
of labels. Probably, methods from geometry or calculus 
would be needed to solve this puzzle (Sec. I VI Bp . (iii) 
In close connection with this, a further geometry-related 
conjecture could be drafted about the nonempty classes: 
They are of nonzero measure. It is known in the bipar- 
tite case that the set of separable states is of nonzero 
measure [H, H3], which can motivate this conjecture, (iv) 
We have given also the necessary and sufficient criteria 
of the classes. This was done by convex roof extension, 
which is a general method for the identification of convex 
subsets, having advantages and disadvantages, (v) First 
of all, convex roof extensions are hard to calculate. How- 
ever, necessary and sufficient criteria for the detection 
of convex subsets seem always to be hard to calculate, 
since they always contain an optimization problem, such 



as finding a suitable witness, or positive map [37[ , or 
symmetric extension |48h50|| , or local spin measurements 
p, or detection vector [5ll. [52|. or local bases [53j], etc. 
(The latter three are for only necessary but not sufficient 
criteria. For further references, see the reference lists of 
[l9l [20j.) These optimization problems have no solutions 
in a closed form in general cases, (vi) Another disadvan- 
tage of convex roof extensions is that this is a "clearly 
theoretical" method, under which we mean that the full 
tomography of the state is needed, then the criteria are 
applied by computer. The majority of the other known 
criteria share this disadvantage. Exceptions are the cri- 
teria by witnesses [13] and by local spin measurements @ 
(only necessary but not sufficient), where the criteria can 
be used in the laboratory, by the tuning of measurement 
settings. However, the optimization still has to be done 
by the measuring apparatus, (vii) An advantage of the 
convex roof extension is that it works independently of 
the dimensions of the subsystems, so the criteria by that 
work for arbitrary dimensions. However, the numerical 
optimization depends strongly on the rank of the state, 
which can be high if the dimension is high, resulting slow 
convergence, (viii) The greatest advantage of the criteria 
given by convex roof constructions — at least for us — is 
that they have a very transparent structure; they reflect 
clearly the structure of the PS classes by construction 
[see (PI)]. 

Now, we turn to some remarks about the three-qubit 
case, which is although particular but very important, 
(ix) First, note that the FTS approach of three-qubit 
entanglement [11] is coming from the famous Black 
Hole/Qubit Correspondence [Hj]. The FTS approach 
has turned out to be fruitful also in some other fields 
of quantum entanglement (55l . [56| . There are some ad- 
vantages of the FTS approach in the three-qubit case, 
although, as we have seen, criteria for the PS subsets 
can be obtained without the use of that, (x) Since the 
convex roof extensions of polynomials can be known to 
be semialgebraic functions [571 . |58| . it can be useful to 
use LU-invariant homogeneous polynomials for the iden- 
tification of the classes. Then we have polynomials of 
this kind from (f?]) coming from the FTS approach, and 
from (|17p with the Tsallis entropy for q — 2 coming from 
the general constructions. The former ones are of lower 
degree, which may lead to more simple convex roof ex- 
tensions, (xi) Moreover, this holds also for the functions 
g a in the general tripartite case if Raggio's conjecture 
holds (Sec. IV Ap . (xii) A little side result of our work is 
that Raggio's conjecture holds for two-qubit mixed states 
which are, at the most, of rank 2. (xiii) An interesting 
question is as to whether all pure-state indicator func- 
tions can be obtained without products of local entropies, 
but using only linear combinations of them, (xiv) We 
note that there are also recent attempts to study the 
general structure of LU invariant homogeneous polyno- 
mials [59l - l63j . Looking for convex roof extensions in the 
language of LU-invariant polynomials would be an in- 
teresting direction of research, (xv) As a disadvantage 
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of the FTS approach, we have to mention that some of 
the indicator functions coming from the FTS approach 
are not nonincreasing on average (|35j) . namely g a and t 
given in (|4d|) and (|4"c| . [Counterexamples for (|3"5|) can be 
constructed for these functions by direct calculation.] 

Finally, we would like to summarize some arguments 
for the relevance of the extension of the Seevinck-Uffink 
classification, (xvi) We can get back the classification 
given by Seevinck and Uffink if we simply forget about 
the sets g € V b \ ac ^ ab , and the functions g^{g). How- 
ever, the appearance of the g a (V>) polynomials is natural 
in the light of the formulas (jtS|) . (j4c]) . and ([id| . This 
motivates the introduction of the sets g € 2?b|ac,c|ab ^ Q 
the classification, (xvii) The g a functions are interest- 
ing in themselves (see Appendix IA 3[) . For all non-GHZ 
vectors, they coincide with the Wootters concurrence- 
squared of two-qubit subsystems (|A10I) . However, note 
that the Wootters concurrences of two-qubit subsystems 
are not suitable for being indicator functions, since they 
can be zero also for GHZ-type vectors, so they do not ful- 
fill the last row of g a columns of Table IIII1 For example, 
for the usual GHZ state (flcj) . the Wootters concurrences 
of two-qubit subsystems are zero, (xviii) In Sec. HV] we 
have shown states which are definitely in classes that are 
different in the extended classification. This is another 
reason for using also the sets g € j)b\ac,c\ab m ^ e c i ass j_ 
fication. 

In closing, there is an important question, which can 
be of research interest as well, (xix) The PS classification 
is about the following issue: "From which kinds of pure 
entangled states can a given state be mixed?" Another 
issue, which is equivalently important from the point of 
view of quantum computation but which we have not 
dealt with, is "Which kinds of pure entangled states can 
be distilled out from a given state?" What can be said 
about the latter? 
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Appendix A: For the three-qubit invariants 



three-qubit state vectors: 

7oW=tr(7r) = |H| 2 , 
l ty)=tr(0, 

h(ip) = 3tr[(7r b (g) 7T c )7r bc 

h{iP) = |Dct(^)| 2 , 



(Ala) 
(Alb) 

tr(7r b 3 )-tr(7r c 3 ), (Ale) 
(Aid) 



where 7r = |^)(-0|. Here I4 is the Kempe invariant [26[ 
(the same for all different b, c £ {1, 2, 3} labels), arising in 
connection with hidden nonlocality. We can alternatively 
form the set of LU-invariant polynomials of (H| as 



n 


= h, 




V 


= 2^o - 


> + 


Sa 


=m 


-la), 


9a 


= I 2 o + 


Ia-h- 


t 


3 


10 t3 


T 2 


= 4/ 5 . 





h + h) 



(A2a) 
(A2b) 

(A2c) 
(A2d) 

(A2e) 

(A2f) 



Obviously, these are not independent. 



2. For the LU canonical form 

In [23|, [24|, the following LU- canonical form ("general- 
ized Schmidt decomposition" ) is obtained for normalized 
pure three-qubit state vectors, 



ml 100) 
%|no) 



fte|101> 

Tfc|lll) 



(A3a) 



with the phase < a < n and the amplitudes T]i > 0, the 
latter being normalized: 



4 

i=0 



(A3b) 



The LU-invariant homogeneous polynomials of (@| for 
this canonical form can be written in a more convenient 
way using another set of independent LU-invariant ho- 
mogeneous polynomials [23L Izjj . which is equivalent to 
that of (|A1[) . This set calculated for the canonical form 
H) is 



In this appendix, we list some features of the LU- 
invariant homogeneous polynomials given in Q for three- 
qubit vectors. 



In the standard basis 



J\ = A, 
J2 = i]QV2, 
J-3 = 110V3, 
•h = 7?0??4, 

^5 = % (A + 772773 



7?17?4l 



(A4a) 
(A4b) 
(A4c) 
(A4d) 
(A4e) 



In [25], the following set of algebraically indepen- 
dent LU-invariant homogeneous polynomials is given for 



where A = | y/rji n^e 10 — ^772773 1 2 . Calculating the polyno- 
mials of (01 for the canonical form (|A3I) . we can identify 
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these invariants as 



y = 4J 4 + -(Ji + J 2 + J 3 ), 


(A5a) 


S a = 4( J 4 + J1+J2+ J 3 ~ Ja) , 


(A5b) 


g a = 2J 4 + 4J a , 


(A5c) 


t = 4J 4 + 8J 5 , 


(A5d) 


r 2 = (4J 4 ) 2 . 


(A5e) 



3. Connections with the Wootters concurrence and 
fidelity of two-qubit subsystems 

The concurrence c(x) for a two-qubit pure state |x) € 
W 1 ®W 2 is given by the concurrence (|12b[) of the one- 
qubit subsystems as 

c(x) = C[tr fe (|x)<xl)] =I<X|X)I, (A6) 

where (x| = -e ® e|x). [Note that e € Lin("H a -> %°*), 

so |x) = (xl = — £ n' £ jj'X l 3 Hj)) leading to the usual 
expression of the spin flip [64] with er 2 = —ie. On 
the other hand, note the familiar expression: (x|x) — 
—Eu'Ejjix 1 3 X % 3 ■} Its square is the local entropy (f2b|) : 
c2 (x) — s a(x) = 4dct(tr b |x)(xl)> this time for the two- 
qubit pure state, a, b € {1, 2}, a ^ b. 

In the convex roof extension of (|A6I) . called Wootters 
concurrence, the minimization can be carried out explic- 
itly [iil [fill , and for a two-qubit mixed state cj it is given 
by the famous formula 

c u (w) = max{0, A j — A2 — A3 — A4}, (A7) 

where A^s are the decreasingly ordered eigenvalues of the 
positive matrix yj y/cuuiy/uj (being the same as the square 
root of the eigenvalues of the non-Hermitian matrix louj), 
and the spin-flipped state is <D = e <g> sui*s ® e. 

If the two-qubit mixed state for which the Wootters 
concurrence-squared is calculated is reduced from a pure 
three-qubit state tt = \ip)(ip\ as uj = ir bc = tr a (|^)(-0|), 
as was investigated in then ir bc is, at the most, of 

rank 2, and [c u (7r &c )] 2 = (Ai-A 2 ) 2 = tr(7r 6c 7r &c ) -2AiA 2 . 
One can check that 

tr(7: bc n bc ) = tr[ 7a (V') t 7a(V')] , (A8) 

which is just g a (ip) [see in (|4d"|) ]. and 

AxA 2 = |det[ 7 .(^)]| = |Det(V)|, (A9) 

(see [Hjl)- The concurrence is then given by 

[ C u (7r, c )] 2 = 5Q (^)-ir(^. (A10) 

The CKW (Coffmann-Kundu-Wootters) equality [3 
about entanglement monogamy, 

Sa W = [c U (7T ab )] 2 + [ C U (^ ac )] 2 +r(^), (All) 



is then equivalent to (l4cl) . 

The roof extension relates the concurrence with an- 
other important quantity, the fidelity (l6j . The fidelity 
between two density matrices u) and a is F(oj, er) = 
tr yj y/ujay/uj, which is the square root of the transition 
probability, and it is in connection with distances and 
distinguishability measures on the space of density ma- 
trices [l|. The fidelity of a state with respect the spin 
flip is F(u),L)), which is just the concave roof extension 
of c(x): 

c n H = F(w,£) = A 1 +A 2 + A 3 + A4 (A12) 

in the two-qubit case. [The concave roof extension is the 
maximization of the weighted average over the decom- 
positions instead of the minimization (|8"|).] Again, for 
the mixed states of two-qubit subsystems arising from 
a three-qubit system being in a pure state, the fidelity 
is [c n (7r bc )] = (Ai + A 2 ) 2 = tr(7r bc 7r hc ) + 2A X A 2 . Using 
(|A8j) and (|A9j) , it is of the form similar to the concurrence 

[c n (7r bc )] 2 = gaW + \rW. (A13) 

Then, using (l4"c]l , we get a CKW- like equality for the 
fidelities, 

Sq (V) = [C^afc)] 2 + [P^ac)} 2 ~ T(i>). (A14) 

On the other hand, from (|A10I) and (|A13I) . <? (^) is just 
the average of the concave and convex roofs and t(i/>) is 
their difference 

9aW = l([c n (n bc )] 2 +[c u (n bc )} 2 ), (A15) 

- [c n kf- [c u k)] 2 . (A16) 

Hence, on the zero-measured set of non-GHZ pure states, 
the convex and concave roof extensions of c(x) are equal, 
and both of them are equal to g a {^)- 



4. The ranges of the invariants 

In this appendix, we show that the LU-invariant poly- 
nomials of Eq. @ range from to 1 if the state vec- 
tor is normalized. [Otherwise, the maximum of the 
functions goes by the corresponding power of the norm: 
< y(i/j),s a {i>),g a (i/j) < n 2 {4,), < t(ip) < n 3 0/>), 
< r 2 (^) < n 4 (^).} 

It is well known that s a (ip) < 1, and this upper 
bound can be achieved. Indeed, using the inequality of 
arithmetic and geometric means (AM-GM inequality), 

a a (1>) = 4det(^ a ) = 4AxA 2 < 4[|(A X + A 2 )] 2 = 1 (for 
the eigenvalues of 7r a ), which is saturated if and only if 
A 4 = A 2 = |, for example, for |B) & ® |0) c , where |B) is 
a Bell state of (TTa|) . 
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From (|4c|) it is clear that g a (ip) < s&(V0 < 1. To show 
that the g a (ip) < 1 upper bound can be achieved, we note 
that, for example, g a (|0)a ® |B)bc) = 1- 

From (|4bj) and (|4c| . one can find that y{ip) = 
l/3(si(V') + SaO/O + s 3 (ip)) < 1. To show that the 
2/(V0 < 1 upper bound can be achieved, we note that, 
for example, for the usual normalized GHZ state (fTcj) 
s a (GHZ) = 1. 

It is known [l2| that t 2 (ip) < 1, and this upper bound 
is achieved; for example, r 2 (GHZ) = 1. 

The proof of the t(ip) < 1 inequality is lengthy but 
straightforward, t is LU invariant, so it is enough to 
calculate (|A5d[) given for the LU-canonical form (| A3|) : 



t(ipsch) = 4?7o[?74 + 4^/772% {VV2V3 ~ V^i^cosa)]. 

First we note that for state vectors of genuine tripartite 
entanglement t(tp) 7^ (see in Table UTI]) . so 770 7^ 0. On 
the other hand, t{i/j) — 4770774 [see in (|A5e|) ]. so 774 decides 
to which class the tripartite-entangled state belongs: If 
l^sch) G V w , then 774 = 0, and if |V>Sch) G V GHZ , then 

For |Vsch> G V w , (770 ^ 0, 774 = 0) we have 



*(V>Sch) = 1677o?72?73 < 16^ (770 + 772 + 77 3 ) 3 
= |(1-^<| = «W)<1, 

where we have used the AM-GM inequality, and the 
(IA3b[) normalization of the state. The first two inequali- 
ties allow equalities for 770 = 7/2 = 773, and Ai = 0, respec- 
tively, which results in a state ^g(|000) + 1 101) + |110» 
being LU equivalent to the usual normalized W state 

m- 

For l^sch) G V GHZ , (770 ^ 0, 774 7^ 0) we have to make 
further distinctions. 

If either 772 = or 773 = 0, then 

*W>Sch) = 4?7o?74 < 4^(7/0 + t? 4 ) 2 



(1 _ m _ m _ ^3)2 < ! = t(GHZ) 



by AM-GM inequality and normalization (|A3bp . The 
inequalities can be saturated for r\\ = r\i — ?/3, in which 
case \ip Sch ) = \GHZ). 
If 771 = 0, then 

i(V'Sch) = 4?7o [774 + 4772773] 

< 4?7o [774 + (772 + 77 3 ) 2 ] < 4?7o [774 + ?72 + m] ■. 

where we used again the AM-GM inequality, and that 
x 2 < x for all < x < 1 with the possibility of equality, 
if and only if x = or 1. 772 + 773 7^ 1 since the state is 
normalized and 770 7^ 0. So we have that the inequalities 
are equalities if and only if 772 = 773 = 0, which results 
that t(^sch) < 4770774, which is maximized by the GHZ 
state, as in the previous case. 



If 771 7^ 0, 772 7^ 0, and 773 7^ 0, then 
*Osch) = 4?7o [774 + 4^/772 V3(V 7 l2V3 - V 7 /i'74COS a)] 

< 4770 [774 + 4772773 + 4^772 ?73»7i m)] 

< 4% [ru + (V2 + mf + \{m + % + m + V4) 2 ] 

where the inequalities are equalities for a — n, and 772 = 
V3 = Vi = V4- This resu lts in t{ip Sc b.) < 4r7 77 4 [l + 8774]. 
Using the normalization (|A3b[) in this case: 770 + 4774 = 1, 
we can find the maximum at 770 = g (5 — v7j . This is 
achieved for the state 



|Vm> = - %/7)|000> + 

^|(l + v/7)i(-|100) + |101) + |110) + |111)) 

For this state t{rp m ) = ±{10 + 1y/7) = 0.52815 • • • < 1 
meaning that this is only a local maximum. 

Appendix B: For the PS classification in general 

In this appendix, we prove some statements used in 
the general PS classification. 

1. Inclusion 

Here we prove (|20[) in the following steps: 
V s3 C V a 



M 



Conv (J (J V 6 C Conv |J (J V 1 

P£P8^P aeaj^a 

UU^U l> 7 

P£P&<P aea^<a 

\J(3 G /3,V<5 ^,3«ea:Ha 
V/3 G (3,3a e a : ft < a 



The equivalences (i) and (v) are by definition (|18bj) 
and p^|) . respectively. 

Equivalence (ii) is the only one where it comes into 
the picture that the story is about quantum states. The 

•1= implication holds, since it is true, in general, that 

ConvB C Conv A <= B C A. However, to the =S- im- 
plication we have to use some special properties coming 
from geometry. Obviously, for the extremal points, 

Extr Conv (J (J V s C Conv (J (J V s 

peps^p Peps^p 

C Conv \J \JV 7 , 
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so 7r e Extr Conv U/3e/3 Us^p ^ 1S a ^ so the element of 
ConvlJ QeQ ,lJ 7 _ <Q 'P 7 . However, n is a projector of rank 
1, so it is extremal also in Conv\J aea U 7 -< Q This 
holds for all 7T, so we have that 



Extr Conv (J \J V s C Extr Conv (J (J V 1 . 

Any A sets of projectors of rank 1 have the property that 
A = Extr Conv A; that is, they are all extreme points of 

their convex hulls, which leads to the =M> implication. 

Equivalence (Hi) comes from set algebra. To see the 

=4- implication, we note that the V" sets are disjoint, so 
every V s on the left-hand side of the inclusion appears on 
the right-hand side as a "P 7 , which means that V/3 € (3, 

V<5 -< /3, 3a G a so that 5 ^ a. To see the impli- 
cation, from the condition 3a so that V s C l^^a^ 7 ' 
but for different 5 ^ /3s there may exist different as, so 
U^/j C U aea U 7r < Q ?> 7 , which holds for all (3 e (3. 

Equivalence (iv) comes from the properties of partial 
ordering: ^ is reflexive on partitions, that is, (3 ^ /3, so 

(iv) 



the =4- implication follows from the 5 — (3 choice. On 
the other hand, ^ is transitive on partitions, which is 

(iv) 

just the <== implication: For all o, if 6 -< p and p -< a 
then S -< a. 



2. Partial order on proper labels 

We show that the relation is a partial order on the 
set of proper labels (|2ip . 

Reflexivity on labels: We need that a -< a, which 
means by definition (TIT))) Va € a 3a' 6 a for which 
a < a 1 . This holds with the a' — a choice, since ^ is 
reflexive on partitions, that is, a -< a. 

Transitivity on labels: Suppose that (3 ^ a. and 7^/3, 
so by definition (fT^) V7 e 7 3/3 6 /3, for which 7^/3, 
and for this /3 3a G ct, for which /3 ^ a. Since ^ is 
transitive on partitions, we have that V7 € 7 3a € a for 
which 7 ^ a, which is 7 X a by definition (|19p 

Antisymmetry on proper labels: Let a and /3 be proper 
labels. Suppose that f3 ^ a and a ^ /3, so by definition 
([HI]) V(3 £ /3 1 e a for which /3 ^ a, and for this 
a 3/3' € j3 for which a ^ /?'. Since -< is transitive on 
partitions, we have that f3 < (3' . This can be true only if 
(3 — (3' , since (3 is a proper label, so we have that (3 < a 
and a ^ /3. Since ^ is antisymmetric on partitions, we 
have that V/3 £E /3, 3a G a for which a = /3, which means 
that (3 C ct. Interchanging the roles of a and /3, we have 
that a C /3. Since C is antisymmetric on sets, we have 
that (3 = a. 



3. For the injectivity of labeling with proper labels 



Let a and (3 be proper labels. Here we prove 
the following steps: 



4=k 



C V a and 2? a C V p 
(3 -< a and a. < (3 
(3 = a. 



Equivalence (i) is the antisymmetry of C on sets, equiv- 
alence (ii) is (|20|) on labels, equivalence (Hi) is the anti- 
symmetry of ^ on proper labels. 



4. For the surjectivity of labeling with proper 
labels 



Here we prove (|23|) in the following steps: 



pa7 _ -pa 



D Q7 C V a and 2?" C V a ~< 
cry ^ ct and a ^ 0:7 
7 ^ ex- 



Equivalence (i) is the antisymmetry of C on sets, equiva- 
lence (ii) is pp)) on labels, (a X 0:7 holds always) equiva- 
lence (Hi) is from the observation that (3 < a and < ol 
if and only if /3/3' ^ a, which can be easily seen from the 
definition (Til 



5. Algorithm generating PS sets 

To obtain an efficient algorithm generating the proper 
labels of all PS subsets, it is necessary to consider the la- 
bels as l-tuples of partitions instead of sets of partitions. 
In this case, a = cti, 0.2, ■ ■ ■ , <%i, so the order of the ele- 
ments is considered to be fixed when an Z-tuple is given, 
and the Oj-s are different for different js. [The (...) I- 
tuple brackets are omitted. Note that, contrary to [f|, 
the lower index of the partitions aj here does not refer 
to the k number of L r sets in aj.] Using ordered struc- 
ture has further advantages beyond the obvious one that 
a computer stores data sequentially, so implementing sets 
would mean additional difficulty. Now the algorithm is 
the following. 

1. [Initialization] Fix an order of the partitions, this 
will define a lexicographical ordering for ^-tuples of 
partitions. This is denoted by <. (This is to avoid 
obtaining an Z-tuple more than once and to make 
the algorithm more optimized.) 

2. [Level 1] Using this ordering, we have all the 1- 
tuples of partitions ordered lexicographically. 
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3. [Induction step: obtaining the I + 1-tuples of par- 
titions (level I + 1) from the Z-tuples of partitions 
(level I)] To every a = oti, 011, . ■ ■ , ol\ Z-tuples (com- 
ing in lexicographically ordered sequence) we have 
to append any such partition aj+i (coming in lexi- 
cographically ordered sequence) that 

(i) ai + i -f, aj and ay a;+i for all j = 1, 2, . . . , I, 
and 

(ii) ai+\ > ai. [Because of the lexicographical or- 
der <, it is enough to consider only the last (Ith) 
partition.] 

Then the resulting a — oti, a^, ■ ■ • , oti, I + 1- 
tuples, and also the partitions in every such I + 1- 
tuple are ordered lexicographically. 

The algorithm stops when no new partition can be ap- 
pended to any of the i-tuples, which comes in finite steps, 
since the number of all the partitions is finite. 



6. Geometric means 

Here we prove that if the functions fij : V — > K 
(j = 1, 2, . . . , q) are non- negative and nonincreasing on 
average, 

Mj(V)>0, (Bla) 

m 

^2pi(J,j(ni) < fij(ir), (Bib) 
i=i 

then their geometric mean 

fl = (HlH2 ■ ■■Vq) 1/q 

is also non-negative (trivially) and nonincreasing on av- 
erage 



M (tt) > 0, 

m 



(We use this for functions defined on pure states, al- 
though the following proof does not use that, so the state- 
ment holds also for functions defined on all states.) 

To obtain this, we will need a Cauchy-Bunyakowski- 
Schwarz-like inequality, for non-negative vectors x^' £ 
W n , x ( p > 0: 

m 

^xfM 2) ---^ ) <l|x (1) IM|x( 2 )|| 9 ...||x^|| 9 , (B3) 

i=l 



where the usual q-norm is 



m 

MI*=[X>? 



1/9 



(B4) 



Indeed, if x^ = for some j, then the inequality holds 
trivially, or else 



m (1) (2) 
llv(l)IL IM 



Si) 



- ||x(l)|| g ||x(2)||/--||xW|| 9 



< 



E 



(Q)\q 



||x(i)||«||xP)||S"'||xW||« 



&)\9 



I|XW||? ||xW||g 



1/'/ 



llx(9)H» 



= 1, 



where the inequality follows from the inequality of the 
arithmetic and geometric means, applied to all terms in 
the sum. 
Using this, 



'Y^Pi^i) = EMM^M 71 *) • ••/»g( 7r i)] 



1/9 



1/9 r i \i 1 /9 

1/9 



1/9 



(B2a) < [mi(tt)] 1/9 [ M2 (tt)] 1/9 . . . 1/9 = 

(B2b) where the first inequality is (|B3j} for sc^ = Jfiffi (^i)] 1/9 
and the second inequality is the condition (IBlbl) . 



I 1/9 
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